Spatial separation in a thermal mixture of ultracold 174 Yb and 87 Rb atoms 
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We report on the observation of unusually strong interactions in a thermal mixture of ultracold 
atoms which cause a significant modification of the spatial distribution. A mixture of 87 Rb and 
174 Yb with a temperature of a few /iK is prepared in a hybrid trap consisting of a bichromatic 
optical potential superimposed on a magnetic trap. For suitable trap parameters and temperatures, 
a spatial separation of the two species is observed. We infer that the separation is driven by a large 
interaction strength between 174 Yb and 87 Rb accompanied by a large three-body recombination rate. 
Based on this assumption we have developed a diffusion model which reproduces our observations. 
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In collisions of ultracold atoms, the strength of the in- 
terparticle interaction is determined by subtle details of 
the interatomic potential [1 j. It is impossible in prac- 
tice to predict the behavior of an ultracold sample ab 
initio and the experimental investigation of a single- or 
multi-species ensemble of ultracold atoms may be full 
of surprises. A prominent example is 133 Cs where after 
years of experimental studies it became clear that with- 
out manipulation using external magnetic fields [2j [3] the 
scattering length is extremely large due to a zero-energy 
resonance [4]. 

A fundamental question related to the interaction be- 
tween atomic species is whether different components 
may occupy the same regions in space or separate spa- 
tially. In trapped quantum degenerate gases this ques- 
tion of miscibility is determined by the interplay between 
single-species and interspecies interactions [5]. For such 
systems, phase separation has been observed in a dual 
species Bose-Einstein condensate with 87 Rb and 85 Rb [6] 
as well as in a strongly interacting Bose-Fermi mixture 
of Li2-molecules and 6 Li atoms [7]. 

In a dual-component gas of thermal atoms, the situa- 
tion is somewhat different because the thermal motion, 
which counteracts interaction-driven spatial separation, 
has to be taken into account. However, for collisionally 
dense samples with a mean-free path much smaller than 
the sample size, motion is diffusive and together with 
a loss process with a spatial and/or density dependence 
this may also lead to a spatial separation of components. 
Such a process, for example, contributes to the stratifi- 
cation of chemical elements inside stars [8]. 

In this Letter, we report on the observation of spatial 
separation in an ultracold thermal mixture of co-trapped 
87 Rb and 174 Yb atoms, indicating large interspecies in- 
teractions. The spatial separation is detected as an exclu- 
sion of the Yb atoms from the region of high Rb density 
at temperatures of a few jaK. 

The trapping potential for the two atomic species is 




FIG. 1. (Color online) Sketch of the trapping geometry for 
174 Yb and 87 Rb showing the trap lasers and magnetic field 
coils (not to scale). Coordinate axes x,y, z are indicated. 



composed of a clover-leaf magnetic trap (MT) and a 
bichromatic optical dipole trap (BIODT), as shown in 
Fig. [T] The diamagnetic Yb atoms are only confined by 
the BIODT created by two superimposed light fields with 
wavelengths of 532 nm and 1064 nm. Both fields are red- 
detuned with respect to the dominant transition in Yb 
at 399 nm [9 j. The total potential for Rb, in contrast, is 
created by the MT, and the repulsive (blue-detuned with 
respect to the dominant 780 nm transition) 532 nm light 
field plus the attractive (red-detuned) 1064 nm field [TO] . 
Thus, the optical potential for Yb is always attractive, 
while it can be tuned from repulsive to attractive for Rb 
by changing the powers in the two light fields. The MT 
and the BIODT have their weak axes aligned along the z- 
direction for maximum overlap of the two atomic clouds. 
A related scheme was recently used in Ref. [TT] . 

In a typical experimental sequence, Yb atoms are 
loaded from a magneto-optical trap, operating on the 
% 3 Pi transition at 556 nm [12 , into the BIODT. 
The BIODT light fields have beam waists of iuio64 ~ 
15 /im and w^s2 ~ 16 /im with initial power levels of 
Pi064 = HO mW and P532 = 3.4 W. The relative posi- 
tion of the two beams is actively stabilized. Next, the 
power of the 532 nm beam is ramped over 30 s to a fi- 
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nal level of P532 = 340 mW. After this ramp, the Yb 
atom number has decreased to « 2 x 10 5 and the Yb 
temperature is « 5 /iK. In this configuration, the trap 
frequencies are (1.05±0.02)kHz radially and (9.5±0.2)Hz 
axially and the calculated trap depth is U/k& ~ 60 /iK, 
where &b is the Boltzmann constant. If not stated oth- 
erwise, all errors represent one standard deviation com- 
bined statistical and systematic uncertainty throughout 
this manuscript. 

While the potential for Yb is ramped to its final value, 
an ensemble of 10 T Rb atoms is prepared in a MT in the 
(F = = —1) state at a temperature of « 1.5 /iK. 

The Rb peak density in the MT is a few 10 13 cm -3 . The 
center of the MT is initially located 0.7 mm from the 
BIODT and the Yb cloud is almost unaffected by the Rb 
preparation process, which takes 35 s. 

Subsequently, the MT is moved to the position of the 
BIODT within 500 ms by applying a magnetic bias. At 
this stage, the effect of the BIODT on the Rb is minimal. 
As the two atomic species make contact, we observe rapid 
thermalization of Yb with the colder Rb cloud to the Rb 
temperature. Heating of the Rb cloud is negligible due 
to the large difference in atom numbers. 

The Rb cloud is then compressed by increasing the 
power of the 1064 nm beam within 100 ms, thus making 
the BIODT attractive for Rb. For Pi 06 4 = 140 mW the 
measured trap frequencies for Rb are (800 ± 40) Hz radi- 
ally and (24±l)Hz axially. The Rb temperature increases 
by adiabatic compression to 2 — 3 /iK and the Rb peak 
density rises up to a maximum value of 3 x 10 14 cm -3 . 
We have confirmed experimentally that the two gases are 
in thermal equilibrium after deformation of the potential, 
and define a temperature Tyb = ^Rb = T. The deforma- 
tion changes the total potential depth for Yb by less than 
10%. No signs of quantum degeneracy are observed. 

The experimentally determined Rb cloud shapes are 
consistent with a model calculation based on a Gaus- 
sian beam profile, measured temperatures and trap fre- 
quencies, magnetic moments, and atomic polarizabilities. 
Densities are assumed to follow a Boltzmann distribution 
in the trap potential. Though the qualitative shape is re- 
produced, the modeling of density is limited by the strong 
dependence of the density on uncertainties in the poten- 
tial. The distribution of the Yb when the Rb cloud is not 
present is consistent with a similar model calculation. 

When the compression, and hence the density of the 
Rb cloud is large, we observe a spatial separation of the 
two atomic clouds, with the Yb bracketing the Rb cloud. 
In Fig. |2j simultaneously taken images of both species 
are shown. Additional time-of-flight images of Rb are 
used to determine absolute atom number and tempera- 
ture. The extent to which the two species separate de- 
pends on the number of Rb atoms, the density distribu- 
tion and the temperature. We do not observe a distortion 
of the Rb cloud caused by the presence of the Yb. Any 
such effect is expected to be small because of the small 
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FIG. 2. (Color online) False color in situ images of 87 Rb and 
174 Yb clouds at T « 3 /iK. For 174 Yb we use absorption 
imaging while for 87 Rb we use dark-contrast imaging [15] . 
The 87 Rb peak density increases from (a) 1.4 x 10 cm to 
(b) 7.3 x 10 13 cm" 3 and to (c) 1.4 x 10 14 cm" 3 . In the radial 
dimension the images are resolution limited. 

Yb density nyb ~ 10 13 cm" 3 . The observed demixing 
is accompanied by a rapid loss of Yb atoms, which is 
attributed to inelastic collisions. No significant loss of 
Rb atoms is observed. A 1/e-lifetime of « 200 ms of the 
Yb cloud has been measured for an initial Rb peak den- 
sity n Rb (0) = (3 ± 1.8) x 10 14 cm" 3 . We attribute the 
loss to three-body recombination of Yb+Rb+Rb, which 
is expected to become important when the interspecies 
scattering is strong [T4HT6] . Two-body loss processes do 
not exist in this system due to the lack of structure in 
the spherically symmetric Yb ground state. Although we 
have observed a spatial separation with the MT switched 
both on and off, all results described in this Letter have 
been obtained with the MT switched on. 

We take two approaches to quantify the spatial sep- 
aration. First, we attempt to explain the phenomenon 
in terms of a mean field interaction. Second, we model 
the system via a diffusive model that includes three-body 
recombination. 

The mean-field approach leads us to quantify the spa- 
tial separation of Yb and Rb in terms of an additional 
potential for Yb due to the elastic interaction with Rb 
atoms. Formally, the total potential seen by the Yb 
atoms is taken to be 

U Yh (r) = V(r) + U YhRh n Rb (f), (1) 

where V(r) is the trap potential due to the BIODT 
and the interaction parameter C/ybRb quantifies the in- 
terspecies interaction. Both f/ybRb and nRb can depend 
on temperature. Effects of the interspecies interaction on 
the Rb density distribution are neglected because of the 
comparatively low density of Yb. 

We fit Eq. [l] to the experimental data assuming that 
the density of Yb follows a Boltzmann distribution, 
e -u Y b(r)/k B T _ Bggideg the relative trap center position, 
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FIG. 3. (Color online) (a) Typical axial 174 Yb density dis- 
tribution for T = (2.5 ± 0.3) /iK and n Rb (0) = (3 ± 1.8) x 
10 14 cm -3 . The solid red curve is a fit of the 174 Yb den- 
sity distribution using a mean field model while the dashed 
black curve is from a diffusion model, (b) Fitted value of 
^YbRb^Rb(O) using a mean field model as a function of nRb(O) 
for measurements at T = 3.2 /iK. The error bars represent 
only fit uncertainties. 



^YbRb is the sole fit parameter because the Rb density 
distribution and the temperature are determined exper- 
imentally. A sample fit is shown in Fig. [3](a). Using this 
method, we obtain /7YbRb/^B = (9.7±5.2) xlO _14 /iK cm 3 
from a series of measurements at T = 3.2 /iK, as shown 
in Fig. [3](b). A spatial separation of Yb and Rb was 
observed between 1.5 /iK and 6.5 /iK. Within experimen- 
tal uncertainties, this analysis results in a temperature- 
independent value for /TybRb- 

In a standard mean-field approach C/ybRb is related to 
a characteristic length aybRb by £/ybRb = 27rft 2 aybRb/M, 
[T7| H8] where /i is the reduced mass of 87 Rb + 174 Yb. 
Using this relationship we obtain a value of aybRb = 
35000±19000ao, with ao the Bohr radius. One is tempted 
to take this value as a measurement of the zero-energy 
interspecies scattering length since the system is well be- 
low the p-wave barrier of 60 /iK |9j [19]. However, for 
large scattering lengths the system is well outside of the 
energy independent regime, and such an interpretation is 
invalid. A next interpretation would be that the length 
is the thermal average of an energy dependent s-wave 
scattering length [18] . We use the effective- range theory, 
which provides an accurate representation of energy de- 
pendence at large scattering lengths [20] [21]. Thermal 
averaging for any zero-energy scattering length yields a 
value for aybRb that is much smaller than the measured 
length. Thus, we conclude that taking the length to be 
a scattering length or a thermally averaged scattering 
length is incompatible with theory. Nevertheless, the fit- 
ted values of t/ybRb^Rb(0) shown in Fig. ^b) are a mea- 
sure for the total strength of the interaction and clearly 
demonstrate its dependence on the Rb density riRb(r). 

While mean-field theory does not provide a satisfactory 
explanation, the results do suggest that the interactions 
between Rb and Yb are strong and that the elastic cross 
section is large. Correspondingly, the mean free path 
of a thermal Yb atom in a cloud of Rb will be much 



smaller than the size of the cloud. Thus we can model 
the trapped system using the time-dependent diffusion 
equation including loss due to three-body recombination. 
The diffusion equation is 

n(r,t) = -V • j(r,t) - K 3 n| b (r)n(r, t), (2) 

where n(r, t) is the Yb distribution, K% is the three-body 
loss rate coefficient, nRb(r) is the Rb density, which will 
be assumed to follow a Boltzmann distribution and is 
independent of time, and the Yb current density is 

= _ fc B Tr(r) + r(r) 

m m 

which describes diffusion and drift. Here r(r) is the spa- 
tially dependent mean collision time and m is the 174 Yb 
mass. 

The mean collision time is determined by Yb+Yb and 
Rb+Yb collisions and satisfies r(r) _1 = TybRb(r) -1 + 
TybYb(r)" 1 . Here rybRb(r)" 1 = (^rei^YbRb) n Rb (r), 
where oybRb is the interspecies collision cross section, 
and v re \ is the relative velocity. The angle brackets rep- 
resent a thermal average. Similarly, self-diffusion is gov- 
erned by rybYb(r)" 1 = (Vrel.YbCYbYb) ™(r). 

By looking for the solution n(r, t) = e _rt n(r), the dif- 
fusion equation becomes an eigenvalue equation for the 
loss coefficient T. The smallest eigenvalue corresponds 
to the longest lived eigenfunction. This eigenfunction 
is the spatial distribution that would be observed in an 
experiment at long times. For K$ = in a trap that 
is deep compared to the temperature, this solution is a 
Boltzmann distribution, nyb e - y ( r )/ k BT ^ w ^] 1 p = 0. 

Rather than solving the full 3-D equation including 
the self-diffusion, which depends on the Yb density and 
is nonlinear, we choose to reduce the equation to 1-D by 
making a few simplifying assumptions and then integrat- 
ing out the transverse coordinates. The reduction is justi- 
fied by the cylindrical symmetry of the cigar-shaped trap- 
ping potential, taken to be harmonic with no axial de- 
pendence in the transverse direction, and the larger radial 
extent of the Rb cloud. As one may assume rybYb( r ) ^ 
TYbRb(r), we substitute the real spatial dependence of 
TybYb( r ) by a functional form that turns r(r) into a sep- 
arable function of radial and axial coordinates and re- 
moves the nonlinearity. With the known Yb+Yb cross 
section aybYb = 7.7 x 10 -12 cm 2 [22], which is constant 
over the range of energies here, rybYb = 5.2 ms for an Yb 
density of 10 13 cm -3 . After the trap geometry has been 
set, the adjustable inputs to the theory are (^rei^YbRb)? 
i\3, and rybYb( r )- Because the interactions appear to be 
stronger than those seen in any other ultracold thermal 
gas, we initially perform a model calculation of diffusion 
using the unitarity limit, aybRb = ^ttH 2 / \jiv re i) 2 . In this 
limit, TybRb can be computed analytically. For values of 
T = 2.5 /iK and n^(0) = 3 x 10 14 cm -3 , corresponding 
to the experimental parameters of Fig. |3](a), one obtains 
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FIG. 4. (Color online) Model calculation of the asymptotic 
lifetime of trapped 174 Yb versus the Rb+Rb+Yb three-body 
loss rate coefficient, K3, assuming a unitarity- limited scatter- 
ing cross section. Inset shows the axial profile of 174 Yb for 
varying K3 (thin blue). Boltzmann distributions for 174 Yb 
(thick red) and 87 Rb (dashed pink) are given. Lines connect 
key points on the lifetime plot to the corresponding distribu- 
tion in the inset. 



TybRb(O) = 5.3 /is. The unitarity-limited value for K3 
is proportional to T~ 2 and at 2.5 /iK it has a value of 
8 x 1(T 25 cm^- 1 [23]. 

Figure [4] shows 1-D solutions in a harmonic trap, giving 
the lifetime 1/r and profile of the longest lived eigenfunc- 
tion for varying K%. For small the Yb distribution is 
Gaussian and V oc K%. With increasing i\3, a kink in the 
lifetime occurs where the Yb is excluded from the Rb. 
The reduced overlap of the clouds results in a weaker 
dependence of lifetime on K%. The next longest-lived 
symmetric diffusion mode has a lifetime that is nearly 
ten times shorter. 

Though it contains significant simplifications, the 1-D 
diffusion model reproduces the experimental density pro- 
file as shown in the dashed black curve Fig.[3](a). The 
curve is calculated using the model of the axial poten- 
tial that is used in the mean field calculations and the 
observed T and riRb(0). The potential has weaker con- 
finement away from the center than the harmonic ap- 
proximation. Both (^rei^YbRb) and K% are used to fit 
the shape and the lifetime of the experimental data. A 
good fit is obtained for K3 = 1.1 x 10 _26 cm 6 s _1 and 
(^rei^YbRb) one fourth of the unitarity limit, correspond- 
ing to a zero-energy scattering length of |aybRb| ~ 500 ao 
within effective range theory. A range of (^rei^YbRb) and 
K3 give similar profiles with lifetimes consistent with ex- 
periment. A more detailed experimental investigation of 
the Rb density-dependence of the lifetime together with 
full 3-D calculations are required to deduce definite val- 
ues for the collisional parameters, which goes beyond the 
scope of the present manuscript. 

As a conclusion, we have observed spatial separation 
in a thermal mixture of 87 Rb and 174 Yb driven by large 



interspecies interactions. A mean-field model cannot ex- 
plain the observations quantitatively, although it sug- 
gests large interspecies interaction. A diffusive model 
including three-body recombination can reproduce the 
experimental observations with physically reasonable pa- 
rameters. To resolve the exact physical origin of the 
spatial separation, further experimental and theoreti- 
cal investigations in a simplified trapping geometry are 
planned. The evidence of strong interactions in this sys- 
tem together with measurements of thermalization rates 
for other Yb isotopes has already been used to predict 
collisional properties and binding energies near thresh- 
old for all isotopic combinations of Yb and Rb [24] . 
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